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The accurate prediction of jet noise is possible only if the jet flow field can be predicted accurately. 

Predictions for the mean velocity and turbulence quantities in the jet flow field are typically the product of 
a Reynolds-averaged Navier-Stokes solver coupled with a turbulence model. To evaluate the effectiveness 
of solvers and turbulence models in predicting those quantities most important to jet noise prediction, 
two CFD codes and several turbulence models were applied to a jet configuration over a range of jet 
temperatures for which experimental data is available. 


Introduction 

The importance of reducing jet noise in both com- 
mercial and military aircraft applications has made jet 
acoustics a significant area of research. 1 The prediction 
of jet noise has two primary elements: the prediction of 
the jet flow field, including turbulence quantities, and 
the prediction of the acoustic field generated by that 
turbulent flow field. The objective of the present work 
is to evaluate the ability of two modem Navier-Stokes 
flow solvers, coupled with several modem turbulence 
models, to predict the mean flow and turbulence prop- 
erties of a supersonic jet plume. This evaluation is 
accomplished by comparison with experimental jet data, 
permitting the weaknesses of the flow solvers and models 
to be identified. The ultimate goal of this research is 
to develop a reliable flow solver applicable to the low- 
noise, propulsion-efficient, nozzle exhaust systems being 
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developed in NASA focused programs. These programs 
address a broad range of complex nozzle geometries op- 
erating in high temperature, compressible, flows. 

The development of a reliable flow solver to be ap- 
plied to noise prediction entails the proper treatment 
of such issues as gas compressibility, variation of the 
turbulence dissipation rate due to heat transfer and three- 
dimensionality in turbulent shear flow (z. e. y turbulence 
anisotropy.) In a round-jet application of the MGB 
acoustics code 2 for jet-noise prediction, Khavaran 3 used 
the PARC aerodynamic code with the two-equation k-e 
turbulence model and standard Launder coefficients 4 to 
predict the jet-plume turbulence structure. The acoustic 
predictions by this Lighthill acoustic-analogy 5 approach, 
with the necessary turbulent length and time scales pro- 
duced by the flow solver, were in reasonable agreement 
with experimental measurements reported by Yamamoto 
et al 6 However, the extension of this noise-prediction 
approach to a wide class of nozzle geometries and wide 
operating ranges of jet total pressure and temperature 
requires the adequate treatment of turbulent compress- 
ibility and anisotropy effects. Such a treatment is the 
goal of the present research. 

Seiner, Ponton, Jansen and Langen 7 previously 
discussed the jet configuration examined here. This 
convergent-divergent nozzle with an exit diameter of 3.6 
in. was designed for an exhaust Mach number of 2.0 
and a total temperature of 1680° F. The acoustic and 
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Fig. 1 Computational domain 


aerodynamic data reported by Seiner et al? covered a 
range of jet total temperatures from 104° F to 2200° 
F at the fully -expanded nozzle pressure ratio. The 
aerodynamic data included centerline and radial mean 
velocity and total temperature profiles. 

Computations were performed independently with 
two CFD codes, ISAAC 8 and PAB3D. 9 Turbulence 
models employed include the k - e model, 10 the 
Gatski-Speziale algebraic-stress model, 11 the Girimaji 
model 12 and each of these models with the Sarkar 
compressibility correction. 13 Centerline and radial 
values of such quantities as mean velocity, Mach 
number and mean temperature are compared with the 
experimental data. 

Overall, the qualitative behavior of the experimental 
data was captured fairly well by both codes and most 
models. For the mean quantities compared here, the 
algebraic-stress models did not always represent an im- 
provement over the k — e model and the compressibility 
correction did not always represent an improvement over 
the uncorrected models. 

Codes and Models 

Both CFD codes used in this investigation were devel- 
oped at NASA Langley Research Center and are intended 
for use in the solution of turbulent, viscous, flows (em- 
ploying turbulence models) about general geometries. 

ISAAC (Integrated Solution Algorithm for Arbitrary 
Configurations) 8 is a finite- volume code of second-order 
accuracy which solves the full Navier-Stokes equations. 
An upwind scheme based on Roe’s flux-splitting scheme 
is applied to the convective terms, diffusion terms are 
centrally differenced and iterations are performed us- 
ing an implicit, spatially split, approximate-factorization 
scheme. 


The PAB3D code 9 solves the Reynolds-averaged 
Navier-Stokes equations with some stream-wise 
derivatives of the viscous terms neglected. Diffusion 
terms are discretized by a central-difference scheme. 
The implicit iteration operator is constructed with the 
van Leer scheme and the explicit terms (such as the 
convective terms) are evaluated with the Roe upwind 
scheme. 

The grid used in these simulations is constructed of 
five blocks, one in the nozzle (61 x 61 mesh points), one 
external to the nozzle (61 x 61 mesh points) and three 
blocks downstream of the nozzle exit (65 x 121, 97 x 121 
and 97 x 121 mesh points). The computational domain 
is shown in Figure 1 . Velocity, pressure and temperature 
boundary conditions consistent with experimental con- 
ditions are imposed at the upstream boundary inside the 
nozzle and a modest ffee-stream flow of M = 0.05 is im- 
posed at all non-wall boundaries except the downstream 
boundary, where a subsonic outflow boundary condition 
is imposed. No-slip boundary conditions are imposed at 
all walls. The flow is assumed to be axisymmetric. 

The turbulence models employed in this investigation 
include both k — e and algebraic-stress two-equation 
models. These models span a range of models usable in 
practical computations today, from the more well-tested 
models to newer models whose capabilities are not yet 
so fully explored. Tests of these models give an accurate 
picture of the current status of flow-field computations 
for jet noise prediction. 

In all computations, the Favre-averaged compressible 
Navier-Stokes equations 8 for the averaged density p and 
the velocity components 

P,t + {puk),k = o (1) 

(pui),t + (puiUj),j = aij,j - G pT i:f )j , (2) 

are employed, where denotes the Favre-averaged 
viscous stress and commas denote differentiation. The 
molecular viscosity p is determined by Sutherland’s law. 

The Reynolds stress, r*j, is modeled by the expres- 

14 

sion J 

JfTii = \pk8 ii -2pC;fJj-[{S ij -\s kk S ij ) 

+ ai -(S ik W kj + S jk W ki ) 

Ss 

- Qs — ( SikSkj — -SkiSkiSij)], (3) 

£$ o 

Sij — (1/2) (ti*j *f and W{j — — 

Ujj) being the symmetric and anti-symmetric parts of 
the mean-velocity-gradient tensor. The coefficients C*, 
04 and 05 are determined by the individual turbulence 
models discussed below. 

The turbulent kinetic energy k and the solenoidal part 
of the dissipation rate e s (the total dissipation e ~ £ s + £ c 
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is the sum of the solenoidal and curl-free dissipations; in 
the case of no compressibility correction, e c = 0) are 
determined from the conservation equations 

(pk),t + {pujk)j = pP-pe 



£ £ £ 
= pC e ifi-^P - pC e 2h~Y 

+ ( (t* + pc;u£-)kj) + *.. ( 5 ) 

\ a £ e S / J 

P = —TijUij is the turbulence production. 

The k — £ model, the simplest of the models discussed 
here, has <24 = 0, a$ = 0 and C* a constant (commonly 
denoted C^) that takes the value 0.096 in the ISAAC 
model and 0.09 in the PAB3D model. In the k and e s 
equations in ISAAC, <r* = 1.0, a £ = 1.3, C c 1 = 1.44 
and C € 2 = 1*83. In PAB3D, the different value C € 2 = 
1.93 is employed. In order to permit integration to the 
wall, the damping functions 


and 

«3 = ^(2 -C 4 ) 2 , 9 = Qcj+Cs-l) (11) 

The model for the pressure-strain correlation of Speziale, 
Sarkar and Gatski 15 is employed to yield the values 
Ci = 6.8, C 2 = 0.36, C 3 = 1.25, C 4 = 0.40 and 
C5 = 1.88. No wall damping of the Reynolds stress 
is required, so = 0. In ISAAC, the e equation is 
changed to give proper log-law behavior with the new 
stress relationship; 14 the expressions /1 = 1, / 2 = 
(1 - exp(-2/ + /5.5) 2 , e = £, and Xw - 0 are used in 
(5). PAB3D employs the same k and e equations as are 
used in the k — e model and does not use the sixth-power 
terms in (7) introduced to make the model agree with 
rapid-distortion theory in the large-shear limit. 

The Girimaji algebraic-stress model 12 also employs 
the Reynolds-stress expression of (3), in this case with 
the coefficients 

c; = -\g u (i 2 ) 

and 





(6) 


a 4 = - 


LI 


rjiL\G\ 


^5 = - 


2L 3 


* ViLlGi 


(13) 


and /1 = l,/ 2 = 1 are employed in ISAAC, where 
y+ = pyu T /fjL ( y is normal to the wall and u r is 
the wall shear velocity) and Rt — ~pk 2 jp£ s . The 
remaining quantities in the e equation are i = £ s — 
p/p \VVk\ 2 and Xw = exp(-(i? T /40) 2 ] [-2 p£ s i + 
1.5 p{£ s - 2 pk/{py+ 2 )) - 1.5 C el £ 8 P)/k. In PAB3D, 
a low Reynolds number model is used in blocks 1 
and 2 (the blocks adjacent to the nozzle walls) where 
/„ = exp[-3.41(l + Rt)~ 2 ), fi = 1, and f 2 = 
1 — 0.3exp(-J?y). i is defined as in ISAAC and 
Xw = 2pp T /p [Sy(Uijj) 2 + Eij* («?,*)]. In blocks 
3, 4 and 5, a high Reynolds number model is used with 
= fi = h = M = and Xu> = 0. 

The Gatski-Speziale algebraic-stress model 11 employs 
the full Reynolds-stress expression (3) with coefficients 

3(1 + rf) 4- 0 . 2 ( t ? 6 4 - C 6 ) _ 

" 1 3 + 77 2 + 6t7 2 C 2 + 6C 2 + V 6 + C 6 ’ 

a 4 = |(2 — C^), a 5 = g{ 2 — C$), (8) 

that depend on the flow quantities 


Here, 771 = k 2 / £ 2 Si jSij and 


= \C° - 1, L\ = Cl + 2, 



The factor G\ is the (explicit) solution of a cubic alge- 
braic equation; see the original work of Girimaji 12 for 
details of this solution and of the selection of the phys- 
ically relevant root. The PAB3D implementation of this 
model again employs the k — £ model e equation. Other 
than the PAB3D version of this model implemented as 
part of this work (see below), this model has not been 
implemented in ISAAC. 

The Sarkar compressibility correction 13 provides for 
the computation of the compressible part of the turbulent 
dissipation e c in terms of the turbulent Mach number 
M t = y/Tii/a 00 and the solenoidal dissipation £ s : 

£ c = a c M 2 £ s (15) 


tj 2 = a^SijSij, C 2 = as^WijWij (9) 

with 

aiaB i(i -Ca )’ q * = t( 2 - c 3 ) 2 ’ < l0 > 


The constant is here taken to be a c = 0.5. 

Results and Discussion 

Computations were performed at jet total temperatures 
of 104° F and 1550° F. Runs were typically made with 
the grid coarsened by a factor of 2 or 4 in each direction. 
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Fig. 2 104° F case, no compressibility correction. Exper- 

iment, solid line; ISAAC k — e model, dash-dot-dot line; 
ISAAC GS ASM, long-dash line; PAB3D k - s model, dot- 
ted line; PAB3D GS ASM, dash-dot line; PAB3D Girimaji 
ASM, dashed line. 

but runs on the full grid were also performed to confirm 
the grid-independence of the results. All runs were 
continued well past the point where no changes at the 
level of graphical accuracy were seen over an extended 
period of iteration. 

The variation along the centerline of the mean veloc- 
ity, Mach number and temperature as computed by the 
models intrinsic to ISAAC and PAB3D are shown in 
Figures 2-5. It is seen that, in most cases, agreement 
is reasonable for at least some of the runs. Oscillations 
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Fig. 3 104° F case, Sarkar compressibility correction. 

Lines as in Figure 2. 

in the velocity and temperature near the nozzle exit are 
present in both the experimental and numerical results. 
Schlieren data presented by Seiner ex ah 1 show a weak 
train of shock waves in the plume; this is normal for a 
nozzle designed for a single jet total temperature. The 
numerical simulations, performed in a computational 
domain that extends into the nozzle interior, also show 
weak shocks and these are likely to be the cause of 
the oscillations. Except for details of these oscillations, 
all runs gave essentially identical results — in good 
agreement with experiment — within the potential core 
(the region of essentially inviscid flow downstream of 
the nozzle between the shear layers emanating from the 
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nozzle walls.) The prediction of the end of the potential 
core region, on the other hand, is a matter of some del- 
icacy, involving as it does the accurate prediction of the 
viscous and turbulent physical mechanisms of the shear 
layers over the entire distance from the nozzle exit to the 
eventual coalescence of the shear layers at the centerline. 
This prediction, not surprisingly, varies in the results for 
the different codes and models discussed here. The re- 
sults downstream of the coalescence of the shear layers, 
while exhibiting some offset (due, presumably, to mis- 
prediction of the coalescence point), have approximately 
the correct slope in most cases. An interesting additional 
discrepancy between experiment and computation is the 
abrupt departure at the end of the potential core from 
the potential core values in many of the computational 
results, contrasted with the smoother transition seen in 
the experiments. While there is no clear superiority of 
the algebraic-stress models in the results shown here, it 
should be kept in mind that a comparison of the details 
of boundary and shear layers could well show advantages 
of the ASM’s. 

The Sarkar compressibility correction does not repre- 
sent a universal improvement in the simulations reported 
here. In the 104° F jet case, the correction makes the 
PAB3D predictions significantly worse, and the ISAAC 
predictions are slightly high by about as much as they 
were low for the computations without the compressibil- 
ity correction. In the 1550° F jet case, the compressibil- 
ity correction significantly improves the PAB3D results, 
but makes the ISAAC results dramatically worse. 

Despite the generally good agreement between the 
present results and experiment, there are noticeable dif- 
ferences between the results of the two codes and the 
different models used in the codes. For the most part. 
Figures 2-5 show a grouping of the results with ISAAC 
and its different models and of PAB3D and its differ- 
ent models. This fact suggest the differences in results 
may be due more to differences in the codes than to 
differences in the models. In order to further investi- 
gate these differences, the k — e model in PAB3D was 
implemented in ISAAC (except for the PAB3D form for 
Xum which involves derivatives not readily available in 
the ISAAC code), as were the PAB3D implementation of 
the Gatski-Speziale and Girimaji ASM’s. Comparisons 
between the ISAAC implementations of the models, the 
PAB3D implementations of the models and the PAB3D 
models implemented in ISAAC are shown in Figures 
6, 7 and 8 for the 104° F case (with no compressib- 
lity correction.) The results ffom ISAAC with its own 
models and with the PAB3D models were very similar 
for the k - e model and the Gatski-Speziale algebraic- 
stress model, indicating that the differences between the 
ISAAC and PAB3D results are primarily due to the codes 
themselves. The Girimaji ASM had not been previously 





Fig. 4 1550° F case, no compressibility correction. Lines 

as in Figure 2. 

implemented in ISAAC, so Figure 8 contains results only 
from PAB3D and from the PAB3D implementation of 
the Girimaji ASM in ISAAC. Curiously, these results 
are significantly closer than the PAB3D and ISAAC- 
with-PAB3D-model results for the k — € model and the 
Gatski-Speziale algebraic-stress model. An interesting 
effect encountered in the implementation of the PAB3D 
models in ISAAC was the convergence to a completely 
unphysical solution due to the differences in the flow- 
field intialization of the two codes. A different initial- 
ization was required in order to make ISAAC with the 
PAB3D models converge to the physical solution. 

Finally, radial experimental data was available for a 
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Fig. 5 1550* F case, Sarkar compressibility correction. 

Lines as in Figure 2. 

900° F case, and a comparison between this data and 
the computational predictions at several positions down- 
stream of the nozzle exit are shown in Figure 9. The 
predictions of ISAAC and PAB3D are very similar, dif- 
fering primarily near the centerline, where differences 
have already been exhibited in the earlier Figures. The 
predictions are in reasonable agreement with the experi- 
mental results. 

In summary, overall agreement of the computations 
with experiment is good. The two codes each gave fairly 
consistent results with the different turbulence models, 
and the differences between the codes seemed to be 
greater than the differences between the models. Addi- 
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Fig. 6 104° F case, comparison of k-e models. Experi- 

ment, solid line; ISAAC, dash -dot line; PAB3D, dashed line; 
PAB3D model in ISAAC, dotted line. 


tional evidence for this was given by the computations 
with ISAAC using the PAB3D versions of the models, 
which gave results much closer to the ISAAC results with 
its own models than to the PAB3D results. Possible rea- 
sons for these differences between the results of the two 
codes include different handling of viscous fluxes, first 
vs, second order advection in the turbulent equations and 
other differences in the numerical algorithms employed 
in the two cases. These differences are currently under 
investigation, as are improvements to the representation 
of the physics of the problem, such as the incorporation 
of a variable Prandtl number. 
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Fig. 7 104° F case, comparison of Gatski-Speziale ASMs. 
Lines as in Figure 6. 
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